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Abstract 

We obtain the exact energy spectra and corresponding wave functions of the radial Schrodinger 
equation (RSE) for any (n, /) state in the presence of a combination of psudoharmonic, Coulomb 
and linear confining potential terms using an exact analytical iteration method. The interaction 
potential model under consideration is Cornell-modified plus harmonic (CMpH) type which is a 
correction form to the harmonic, Coulomb and linear confining potential terms. It is used to 
investigates the energy of electron in spherical quantum dot and the heavy quarkonia (QQ-onia). 
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I. INTRODUCTION 



The problem of the inverse-power potential, 1/r", has been used on the level of both 
classical and quantum mechanics. Some series of inverse power potentials are applicable 
to the interatomic interaction in molecular physics [1-3]. The interaction in one-electron 
atoms, muonic, hadronic and Rydberg atoms takes into account inverse-power potentials 
[4]. Indeed, it has also been used for the magnetic interaction between spin- 1/2 particles 
with one or more deep wells [5] . The analytical exact solutions of this class of inverse-power 
potentials, V{r) — Ar~^ + Br~^ + Cr~'^ + Dr~^, A > 0, were presented by Barut et al. 
[6] and Ozgehk and §im§ek [7] by making an available ansatz for the eigenf unctions. The 
Laurent series solutions of the Schrodinger equation for power and inverse-power potentials 
with two coupling constants V{r) = Ar^ + Br~^ and three couphng constants V{r) = 
Ar^ + Br^^ + Cr^^ are obtained [8,9]. 

The analytic exact iteration method (AEIM) which demands making a trial ansatz for the 
wave function [7] is general enough to be applicable to a large number of power and inverse- 
power potentials [10]. Recently, this method is applied to a class of power and inverse-power 
confining potentials of three coupling constants and containing harmonic oscillator, linear 
and Coulomb confining terms [11]. This kind of Cornell plus Harmonic (CpH) confining 
potential of the form V{r) = ar'^ -\- br — cr~^ is mostly used to study individual spherical 
quantum dots in semiconductors [12] and heavy quarkonia (QQ-nia) [13,14]. So far, such 
potentials containing quadratic, linear and Coulomb terms have been studied [15,16]. 

The present work considers the the following confining interaction potential consisting of 
a sum of pseudoharmonic, linear and Coulombic potential terms: 

c d 

V{r) = l^har(r) + Vcorn-mod(r) = or^ + 6r - - - — , a > 0, (1) 

where a, b, c and d are arbitrary constant parameters to be determined later. The above 
potential includes the well-known funnel or Cornell potential, i.e., a Coulomb plus Linear 
static potential (CpH), Vcom('") = br — c/r [13], and a term —d/r"^ is incorporated into the 
quarkonium potential for the sake of coherence [14]. We will refer to the potential model 
(1) as a Cornell-modified plus harmonic (CMpH) potential, since the functional form has 
been improved by the additional —d/r"^ piece; besides the contribution from the additional 
term also alters the value of b and c [14,17]. The authors of Refs. [14,18] did not consider 
the harmonic or power-law as the results are expected to be similar. The CMpH potential 
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is plotted in Figure 1 for the values of parameters: a — 1 eV.fm ^, b — 0.217 eV.fm ^, 
c = 0.400 eV.fm and d = 0.010 eV.fm^. 

We will apply the AEIM used in [7,11] to obtain the exact energy eigenvalues and wave 
functions of the radial Schrodinger radial equation (RSE) for the CMpH potential for any 
arbitrary (n, /) state. 

The paper is structured as follows: In Sect. 2, we obtain the exact energy eigenvalues 
and wave functions of the RSE in three-dimensions (3D) for the confining CMpH potential 
model by proposing asuitable form for the wave function. In Sect. 3, we apply our results 
to an electron in spherical quantum dot of InGaAs semiconductor. The relevant conclusions 

are given in Sect. 4. 



II. EXACT SOLUTION OF RSE FOR THE CONFINING POTENTIAL MODEL 



The three-dimensional (3D) Schrodinger equation takes the form [19] 



——A + V(r) 



(2) 



with 



^^2 9_ L\e,(p) 



where m is the isotropic effective mass and Eni is the total energy of the particle. For any 
arbitrary state, the complete wave function, ■0(r, 6, </?), can be written as 



V^(r,^,(^) = ^Ar,^„,(r)l;^(^,^), 
where spherical harmonic Yi^iO, ip) is the eigenfunction of 1/^(6*, ip) satisfying 

L^ie, ip)Yimie, if) = 1(1 + i)n^Yi,nio, v), 

and the radial wave function ip^^i (r) is the solution of the equation 

2d l{l + V 



(3) 



)2m 
^n/(r) + ^[K.-nr)]^n.(r)=0, 



(4) 



(5) 



dr'^ r dr r 

where r stands for the relative radial coordinates. The radial wave function i^nii'"') well- 
behaved at the boundaries (the finiteness of the solution requires that ■0n/(O) — '^nii'"' ~^ 
oo) = 0). Now, the transformation 



(6) 



reduces Eq. (5) to the simple form 



'-'nl 



r + 



en,i - a^r^ - hr + j + 



V) = 0, 



where (l>ni{r) is the reduced radial wave function and 



2m 



2m 



2m, 



2m 



2m 



d. 



(7) 



(8) 



The analytic exact iteration method (AEIM) requires making the following ansatze for the 
wave function [9], 

(Pniir) ^ fn{r)exp[gi{r)], (9) 



with 



fn{r) 



1, n = 0, 

.n(r-af^), n=l,2,- 



gi[r) = --ar^ - I3r + 5\nr, a > 0, (3 > 0. 



(10a) 



(10b) 



It is clear that fnif^) are equivalent to the Laguerre polynomials [20]. Substituting Eq. (9) 
into Eq. (5) we obtain 

/;:(r) + 2ft'(r)/;(r)- 



€iir)=(9'i'{r) + gl'{r) + 



fn{r) 



(11) 



and comparing Eq. (11) and Eq. (7) yields 
Ci l{l + 1) - di 



air + &ir h 

r 



(12) 



First of all, for n = 0, let us take /o(r) and gi{r) given in Eq. (10b) to solve Eq. (12), 



air^+6ir-£or 



Ci /(/ + 1) - di 2 2 



= a^r^+2a^r-a [1 + 2 (5 + 0)]+/3 



2 2/3(5 + 0) 5(5-1) 



(13) 

By comparing the corresponding powers of r on both sides of Eq. (13) we find the following 
corresponding energy and the restrictions on the potential parameters. 



J = ai > 0, 

Ci = 2/3(5 + 0), 



(14a) 

(14b) 
(14c) 
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^ (1 ± /') , where /' = J {21 + 1)' - (14d) 
2 V ^ 

eoi^a[l + 2{5 + 0)] - /3l (14e) 



Actually, to have well-behaved solutions of the radial wave function at boundaries, namely 
the origin and the infinity, we need to take 6 from Eq. (14d) as 

5 = 1(1 + 0. (15) 

Therefore, the lowest (ground) state energy from Eq. (14e) together with Eqs. (14a)-(14c), 
Eq. (15) and Eq. (8) is given as follows 

^"^Vl^sP + '^-FITW '''' 

where the parameter c of potential (1) should satisfy the following restriction: 



c = 



(17) 

Furthermore, the substitution of a, /3 and S from Eqs. (14a), (14b) and (15), respectively, 
together with the parameters given in Eq. (8) into Eqs. (9) and (10), we finally obtain the 
following ground state wave function: 



2ma 9 2mc 
- 



r r (1 + /') 



(18) 




exp 



2m mc^ 
2 V ^h^l + l'Y 



where Di,{z) are the parabohc cyhnder functions [21]. It should be noted that the 
above solutions are well-behaved at the boundaries, i.e., a regular solution near the ori- 
gin could be 4>ni{r — )■ 0) — >■ r^^+^'y^ and asymptotically at infinity as </'„;(r — )■ oo) — >■ 
exp (— ar^ — /3r) — )■ 0. When 6 = (c = 0), the problem turns to become the com- 
moly known pscudoharmonic oscillator (p.h.o.) interaction (a = muj'^/2), and consequently 
a — mu, (3 = b/co and c = (bS/muj) yielding Eqi — {2 + I') ^ — >2^r^m2 and wave function 



Mr) = A^^r(-i+'')/2exp (-^^r^ - ^^r) , where 



.\ ^74 / 2 



The formula (17) is a relationship between parameters of the potential a, b, c and d. There- 
fore, the solutions (16) and (18) are vahd for the potential parameters satisfying the re- 
striction (17). Moreover, the relation between the potential parameters (17) depends on the 
orbital quantum number / which means that the potential has to be different for different 
quantum numbers. In applying the AEIM, the obtained solution for any potential is always 
found to be subjected to certain restrictions on potential parameters as can be traced in 
other works (see, for example, [7-9,11]). 

Secondly, for the first node (n — 1), using /i(r) — (r — a^'^) and gi(r) from Eq. (10b) to 
solve Eq. (12), 

9 , Cl 1(1 + 1) — dl 9 9^^ 

oir^ + hir -Sii - — + ^ = a^r^ + 2al3r 



- a [1 + 2 (5 + 1)] + - 



I3{5 + 1) + aa\ 



(1) 



+ 



(19) 



The relations between the potential parameters and the coefficients a, p, 5 and q;|^^^ are 



6=1(1 + 1'), eu = a[l + 2{5 + l)]-f3'. 



Cl - 2^ (5 + 1) = 2aoL\^, (ci - 2^5) a^' = 25, 



(1) 



where c\ and (x\' are found from the constraint relations. 



2 I Irna 



(2 + + 



(20) 



(21a) 



4a 



62 ^ (1 + I') 



16a 



2ma 



The energy eigenvalue is 



6 = 2 



2ma (2 + V) c 
(1 + 0(3 + 



1 + 



mc^ V 2m 



7r-(l + 0-l 



^1 + (3 + 
(2 + 0' 



and the wave function is 



1 / 2ma 



exp 



-r — 



(21b) 



, (22) 



(23) 
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with 



r(/' + 2)5i + apr(/052-2(^2y^) aS')r(/' + 1)^3 



where 



2a V a 



2a V a 



and a^i ^ is given in Eq. (21b). If there is a p.h.o. interaction, the energy becomes 

Ell = (4 + r) — 



2 2W' 



(24) 



and the wave function 



with 



(^)''''exp(-iS) 



2maj 



2^1 \ / / 2^ 6^1 \ // 2^1 6/1 

5i = D_(i,+2) \\ , -^2 = D_i, \ , -^3 = ^-(i'+l) 



where 



6 = 2 



rmj (2 + /') c 
(1 + V) (3 + V) 



(25) 



and cy^ = 

Following the analytic iteration procedures for the second node {n = 2) with /2(r) — 
(r — aj^ )(r — ) and gi{r) as defined in Eq. (10b), we obtain 

2 , Cl /(/ + 1) - fii 2 2 

air + bir — S2,i 1 ^ = a r + 2apr 



- a [1 + 2 (5 + 2)] + - 



P(5 + 2) + aJ2(^i 



(2) 



+ 



6{6-l) 



(26) 
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The relations between the potential parameters and the coefficients a, /3, 5, a^^^ and ctg^^ are 



c, - 2/3 (5 + 2) = 2a 5] af\ {c, - 2(36) «f ^ = 2^ E 



i<j 



1=1 



[ci -2^(5 + 1)] J2 «f ^ = 4a J] af^af +2(25 + 1), 
The coefficients ct^^ and q;2^^ are found from the constraint relations, 



i=l 
2 



i=l 



(2)2 



V j=l / i<fc j<k 

Hence, the energy eigenvalue is 



E21 — 



and the associated wave function is 



m 



(27) 



(28a) 



(28b) 



(29) 



(30) 



where af^ and a^^ should satisfy the restriction relations (28a) and (28b). 

We apply the present method for the third node (n — 3) by taking fsir) — (r — af^){r 
Q;2^^)(r — ctg^^) and gi{r) as defined in Eq. (10b) to obtain 



air + bir - £3,; h 



- a [1 + 2 (5 + 3)] + - 



{5 + 3) + aY, 



a. 



(3) 



1=1 



+ 



5{5-l) 



(31) 



The relations between the potential parameters and the coefficients a, /3, 5, af\ a^^ and 

(3). 



dig 'are 



5=\{l + l'), e,,i^ a[l + 2{5 + ?,)]- ^\ 
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3 3 3 

ci - 2/3 (5 + 3) = 2a 5^ af\ {c, - 2^6) ^ «f ^af = 25 ^ af^af\ 

i=l i<j<k i<j 

3 3 

[ci - 2^ (5 + 2)] J] = 4a J] af af + 3 (25 + 2) . (32) 

1=1 i<j 

The coefficients a^^\ a^^ and ag^^ are found from the constraint relation, 

E«P + /9E«r-3(5 + l) = 0, (33) 



a 

j=l i=l 



The energy eigenvalue is 



and the wave function is 



V,3,(r) = A',;ff (r - a."") r'-'+''>/^ exp ('-i \/^r^ - . (35) 

We can repeat this iteration procedures several times to write the exact energies of the 
CMpH potential for any n state as 



B„, = ,/||(2 + 2„ + r)-^, (36) 



and the wave functions is 



^niir) - iV.n (r - af)) r(--^'')/^ exp [-^^r^ - ^J^^br^ • (37) 
The relations between the potential parameters and the coefficients a, (3, 6, a^"\ C(2^\ ■ ■ ■ , 



(n) 

ah. are 



a = V^, P = 7^, 5=^(1 + 0, £2,^ = a[l + 2(5 + n)]-/3^ 
ci -2p{5 + n)^ 0, (n = 0) 

n 

ci - 2/3 (5 + n) = 2a E \ n = 1,2,3,- •• 



i=l 



[ci - 2^ (5 + n - 1)] E = n (25 + n - 1) , (n = 1) 

i=l 



[ci - 2/3 (5 + n - 1)] E = 4" E ^ + n (25 + n - 1) , n = 2, 3, 4, 

i=l i<j 



[ci - 2/3 (5 + n - 2)] ^ ) = (n - 1) (25 + n - 2) ^ af \ (n = 2) 

i<j i=l 

n n 

( ci - 2P5) J2 = 25 \ (n = 3) , 

i<j<k i<j 



n n 

j 

i<j i<j 



[ci-2/3{6 + n-2)]J2 "l"^"^ = - 1) (25 + n - 2) ^ af^aj. 



n 



+ 4a ^ n = 3,4,5,-.-, (38) 

i<j<fe 



and so on. 



III. RESULTS AND DISCUSSIONS 

Now, we consider a special case of potential (1) and an application to our results. For 
example, when 6 = then leads to c = 0, then we have the p.h.o potential, i.e., Vph{r) = 
\muP'r'^ — ^, hence, the energy difference between the ground state and the excited states 
is 

= En - E,i = (4 + r) ^ - (2 + /') ^ = hco, (39) 

which can be used to calculate the values of the potential parameters for the desired system. 

We now apply the present results to describe a realistic physical system called indium 
gallium arsenide (InGaAs) quantum dot, i.e., a piece of this material of a spherical form 
which is considered as a semiconductor composed of indium, gallium and arsenic [11]. It 
is used in high-power and high-frequency (say, uj ~ 10^^ Hz) electronics because of its 
superior electron velocity with respect to the more common semiconductors silicon and 
gallium arsenide. InGaAs bandgap also makes it the detector material of choice in optical 
fiber communication at 1300 and 1550 nm. The gallium indium arsenide (GalnAs) is an 
alternative name for InGaAs. In Fig. 2, we plot the ground state electron energy 

£;oKu;) = ( 2 + (2Z + 1)^ - — - ^ 1 + y {21 + lf-—d\ , (40) 



versus u in the interval 2 x 10^^ < a; < 10 x 10^^ Hz taking the value of c = 0.001 eV.fm 
and d — eV.fw? for the cases I — and I — 1, respectively (harmonic. Coulomb and hnear 
combination terms). In Fig. 3, we take instead the value of the parameter d = 0.1 eV.fm'^ 
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(pseudoharmonic, Coulomb and linear combination terms). The effective mass of electron 
in the InGaAs semiconductor has been chosen as m = O.OSmg and k — 6.5821 x 10~^^ eV.s. 
It is seen from Fig. 2 and Fig. 3 how the increase in the value of u leads to an increase in 
the energy of electron. The flexibility in the adjustment of the parameter d allows one to 
fit the spectrum of the desired model properly (cf. Fig. 2 and Fig. 3). The parameter d 
should satisfy the condition d < {21 + 1)^ h^/ (8m). In Fig. 4 we plot the ground state wave 
function ■0o,i('^) CpR potential for the cases I — and I — I, respectively, using the 

values of potential parameter c = 0.001 eV.nm for an electron with effective mass m — 0.05 
me and frequency a; = 10 x 10^^ Hz. Further, in Fig. 5 we plot the ground state wave 
function ipQi{r) of the CMpH potential for the cases / = and / = 1, respectively, using the 
values of potential parameters c = 0.001 eV.nm and d = 0.01 eV.nm'^ for an electron with 
effective mass m = 0.05 rrie and frequency a; = 10 x 10^"^ Hz. In Figs. 6 and 7, we show 
electron energy as a function of parameter c in the interval 6 x 10~^ < c < 10 x 10"^ eV.nm 
and d — 0.01 eV.nm^ for frequency a; = 8 x 10^^ Hz and effective mass m = 0.05 me for the 
cases I — and I — 1, respectively. From Fig. 5, the increase in c leads in the decrease in 
the electron energy in the InGaAs semiconductor. In Fig. 8, we plot the first excited state 
electron energy 

huj 2 (2 + Vf mc^ 



/ (1 + /') ("i + V) f K^u , , , \ 

versus u in the interval 2 x 10^^ < a; < 10 x 10^^ Hz taking the value of c = 0.001 eV.fm 
and d — eV.fw? for the cases I — and I — 1, respectively. In Fig. 9, we take the value of 
the parameter d = ^.\ eV.fm?. We remark that the strongly attractive singular part —d/r'^ 
is physically incorporated into the quarkonium Cornell potential as the first perturbative 
term for the sake of coherence to describe the heavy quarkonia (QQ-nia) (see, for example, 
[13,14] and the references therein). It also resemles the centrifugal barrier term l{l + l)/r^ 
in the Schrodinger equation. This attractive term —djr^ together with the h.o. part ar^ 
constitute the so-called p.h.o. when 6 = (^ = 0) in Eq. (14b) leading to c = in Eq. 
(14c). 

In Table 1, we calculate the lowest in = 0) energy states (/ = 0, 1 and 2) from Eq. (36) 
and from the numerical solution of the radial Schrodinger equation (7) using the values of 



11 



parameters given by Ref. [34] using the supersymmetry quantum mechanics (SUSYQM). It 
is clear that the calculated energy states in the present work are in good agreement with 
the results obtained numerically and SUSYQM [34]. The accracy of our numerical results 
is 0.0070% - 0.0095%. 

IV. CONCLUSIONS AND OUTLOOK 

In this work, we explored the analytical exact solution for the energy eigenvalues and 
their associated wave functions of a particle in the field of Cornell-modified plus harmonic 
confining potential. We have used the analytical exact iteration method (AEIM) which 
required making a trial ansatz for the wave function. The general equation for the energy 
eigenvalues is given by Eq. (36) with some restrictions on the potential parameters. If one 
takes b = then c = 0, hence, the potential (1) turns to the p.h.o. potential with energy 
eigenvalues: 

= (2 + 2- + y^(2/+l)^-^c^j (42) 

and wave functions: 

^niir) = iV,n (r - aS")) r^^^'')/^ exp (^-^-^r^^ . (43) 

The present results in Eqs. (42) and (43) coincide with Eqs. (15) and (16) of Ref. [22] 
obtained by the exact polynomial method, Eqs. (72) and (78) of Ref. [23] obtained by the 
Nikiforov-Uvarov method and Eqs. (30) and (31) of Ref. [24] obtained by the wave function 
ansatz method after setting Dq/tq — a, DqTq — —do and 2Do — 0. The model solved in the 
present work can be used in modeling the quarkonium [14] perturbed by the field of p.h.o. 
or electron confined in spherical quantum dots [11]. Finally, our solution to this confining 
potential is being considered important in many different fields of physics, such as atomic 
and molecular physics [25,26], particle physics [13,27,28], plasma physics and solid-state 
physics [29-33]. 
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TABLE I: Lowest (n = 0) energy spectra (for % = m = 1). 



a b c dl Numerical Present SUSYQM [34] 
5^ 1 4 -7.618 -7.625 -7.625 
^ 1 8 1 -7.368 -7.375 -7.375 
4 1 12 2 -7.120 -7.125 -7.125 



FIG. 1: A plot of the CMpH potential [see Eq. (1)] with the selected values of parameters: a = 1 
eF./m-2, h = 0.217 eV.fm-^, c = 0.400 eV.fm and d = 0.010 eV-fm"^. 



FIG. 3: The ground state electron energy in InGaAs semiconductor versus u in the field of the 
CMpH potential with c = 0.001 eV.nm and d = 0.1 eV.nm? for the cases I = and I = 1, 
respectively. 



FIG. 4: Behaviour of the ground state wave function ^„^o,i=o('") (dashed line) and ^„=o,/=i(^) 
(continuous line) in the field of the CpH potential with the value of c = 0.001 eV.nm for an electron 
with effective mass m = 0.05 rUe and frequency a; = 10 x 10^^ Hz in the InGaAs semiconductor. 



FIG. 2: The ground state electron energy in InGaAs semiconductor versus oj in the field of CpH 
potential with c = 0.001 eV.nm for cases / = and / = 1, respectively. 
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FIG. 5: Behaviour of the ground state wave function ^„^Q;^Q(r) (dashed Hue) and V'ra=o«=i('') 
(continuous hne) of the CMpH potential with the values of c = 0.001 eV.nm and d = 0.01 eV.nm? 
for an electron with an effective mass m = 0.05 rUe and frequency a; = 10 x 10^^ Hz in the InGaAs 
semiconductor. 



FIG. 6: Ground state energy of electron versus c, for the case I = 0, co = 8x 10 Hz and d = 0.01 
eV.nm?. 



FIG. 7: Ground state energy of electron versus c, for the case / = 1, a; = 8 x 10 Hz and d = 0.01 
eV.nm^. 



FIG. 8: The first excited state electron energy in InGaAs semiconductor versus u> in the field of 
CpH potential with c = 0.001 eV.nm for cases I = and 1 = 1, respectively. 



FIG. 9: The first excited state electron energy in InGaAs semiconductor versus w in the field of 
the CMpH potential with c = 0.001 eV.nm and d = 0.1 eV.nvn? for the cases I = and 1 = 1, 
respectively. 
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